Surface state scattering by adatoms on noble metals 
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When surface state electrons scatter at perturbations, such as magnetic or nonmagnetic adatoms 
or clusters on surfaces, an electronic resonance, localized at the adatom site, can develop below 
the bottom of the surface state band for both spin channels. In the case of adatoms, these states 
have been found very recently in scanning tunneling spectroscopy experimental- for the Cu(lll) 
and Ag(lll) surfaces. Motivated by these experiments, we carried out a systematic theoretical 
investigation of the electronic structure of these surface states in the presence of magnetic and 
non-magnetic atoms on Cu(lll). We found that Ca and all 3d adatoms lead to a split-off state at 
the bottom of the surface band which is, however, not seen for the sp elements Ga and Ge. The 
situation is completely reversed if the impurities are embedded in the surface: Ga and Ge are able to 
produce a split-off state whereas the 3d impurities do not. The resonance arises from the s-state of 
the impurities and is explained in terms of strength and interaction nature (attraction or repulsion) 
of the perturbing potential. 



I. INTRODUCTION 

At crystal surfaces the symmetry is lowered: The three 
dimensional periodicity in the bulk is lowered to the two 
dimensional periodicity at the surface. This leads to the 
occurrence of two-dimensional surface states^, which 
are spatially confined to the surface, since their wave 
functions decay rapidly into the crystal and are strongly 
damped in the vacuum. Surface states can exist only 
in regions of the two-dimensional Brillouin zone, where 
bulk Bloch states are not allowed. They are character- 
ized by a two-dimensional Bloch vector k\ \ in this surface 
Brillouin zone, which describes the propagation in the 
surface plane. A projection of the bulk band structure 
to the surface plane can result in fen -regions, where bulk 
states are forbidden. In these gaps of the projected bulk 
band structure surface states can occur provided that 
their energy is lower than the work function. These two 
conditions guarantee that the wave functions of the sur- 
face states decay exponentially both into the crystal and 
into the vacuum region. 

Recently a strong interest arose concerning an impor- 
tant physical effect associated with the interaction of 
a two dimensional surface state with the states of an 
adatom. It was shown that for Cu adatoms on Cu(lll) 
a bound state splits off from the bottom of the Cu(lll) 
surface stated. This effect was basically predicted by 
Simon£ who stated that in two-dimensional free space 
any attractive potential has a bound state. Gauyacq et 
al& suggested that an adatom-induced localization of the 
surface state may be observed in STS as a peak appear- 
ing just below the surface band edge when a Cs adatom 
is deposited on Cu(lll). 

Using low-temperature scanning tunneling spec- 
troscopy (STS), Limot et al? investigated silver and 
cobalt adatoms on Ag(lll) as well as copper and cobalt 
adatoms on Cu(lll). The bound state appears both for 
magnetic and nonmagnetic adatoms. Moreover, using 
a Newns- Anderson model the authors explained the re- 
sults as arising from the coupling of the adatom's or- 



bital (which was supposed to be the s-orbital) with the 
surface-state electrons, and being broadened by the in- 
teraction with bulk electrons of the same energy. On 
the other hand, Olsson et alA used the same type of 
experiment and performed pseudopotential calculations 
for single Cu adatoms on Cu(lll). The calculated lo- 
cal density of states (LDOS) exhibits several adatom- 
induced peaks. Two of them are assigned to resonances 
deriving from the d z 2 atomic orbital and sp z hybrid or- 
bitals. The third one corresponds to a localization of the 
surface state at the adatom without a specified orbital 
origin. In fact, these adatom-induced peaks appearing at 
the bottom of the surface state were already observed by 
an other experiment- , again without a clear assignment 
of their origin. Davis et al& have observed similar lo- 
calized states in STS measurements for Cr atoms in the 
surface layer of Fe(OOl) surface. Ab-initio calculations 
of Papanikolaou et al? confirmed this and showed that 
similar localized states occur for many other impurities 
in the Fe(OOl) surface. 

In this work we study the origin and the condition 
of existence for such an impurity-induced split-off state 
on the Cu(lll) surface. For this purpose we have per- 
formed ab-initio calculations using the Korringa-Kohn- 
Rostoker (KKR) Green function method for impurities 
on the Cu(lll) surface. We consider single impurities 
of the 3d and Asp elements as adatoms on the Cu(lll) 
surface and as impurities in the first layer. We find that 
split-off states can appear both for adatoms on the sur- 
face as well as for substitutional impurities in the surface. 
In the case of magnetic impurities these states always ap- 
pear in both spin channels and show a very small spin 
splitting. 



II. COMPUTATIONAL ASPECTS 

In the present work we use the full-potential KKR 
Green function method 10 . This method is ideal for treat- 
ing systems involving impurities on or in surfaces and 
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in bulk crystals. Within the KKR method the impuri- 
ties are described by considering a cluster of perturbed 
atomic potentials which includes the potentials of the im- 
purities and the perturbed potentials of several neighbor 
shells. Also in the vacuum region the space is filled by 
cellular potentials, of which the ones close to the impu- 
rity are perturbed. The impurity potential and the per- 
turbed potentials of the neighboring cells are embedded 
in an otherwise ideal unperturbed surface. 

The KKR method is based on multiple-scattering the- 
ory. For non-overlapping potentials the following an- 
gular momentum representation of the Green's function 
G(r + R„, r' + R„«; E) can be derived: 

G(r + R„, r' + R„,; E) = -iVEj2 R n L (r < ; £)tf£(r> ; E)S nn , 

L 

LL' 

Here E is the energy and R n , R n / refers to the atomic 
positions. r< and r> denote respectively the shorter and 
longer of the vectors r and r' which define the position in 
the Wigner-Seitz (WS) cell centered around R„ or H' n . 
The R L (r\E) and H^{r;E) are respectively the regular 
and irregular solution of the Schrodinger equation. 

The structural Green functions G™ L , (E) are then ob- 
tained by solving the Dyson equation for each spin direc- 
tion 

G n L i,(E) = Gli,{E) 

+ Gii:,(E)Ati;, L ,4E)Gi;,?i,(E) (2) 

n",L"L"' 

The summation in (J2J is over all lattice sites n" and 
angular momenta L", L'" for which the perturbation 

M('„ L „,{E) = tl„ L ,,,{E) - t r l!, L ,„(E) between the t ma- 
trices of the real and the reference system is significant. 

o / 

GYi' are the structural Green function of the reference 
system, i.e. in our case the ideal Cu(lll) surface. 

We have carried out our calculations in the local spin 
density approximation (LSDA) with the parameters of 
Vosko et al. 12 . Angular momenta up to l m ax = 3 are 
included in the expansion of the wave functions and up 
to 2l max = 6 in the charge density expansion. We have 
checked these cut-offs to be adequate for our purpose. 

First, the surface Green functions were determined for 
the surface of Cu(lll). The lattice LSDA equilibrium 
parameter was used (6.63 a.u. « 3.51 A). To describe 
the impurities on the surface, we consider a cluster of 
perturbed potentials which includes the potentials of the 
impurities and the perturbed potentials of several neigh- 
bors shell, with typical size of 29 perturbed sites for the 
adatoms. All adatoms are assumed to sit at the hollow 
position in the first vacuum layer. In the following, we 
take as reference energy the Fermi level Ep. 



According to the Tersoff-Hamann modeli^ the scan- 
ning tunneling spectra can be related to the s-DOS in- 
duced by the surface or by the adatom at the position 
of the STS tip. Adopting this model, we calculate the 
s-LDOS at a distance z — 2.86 A directly above the 
adatom. This corresponds to a lattice position in the 
third vacuum layer above the surface. 



III. RESULTS 
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FIG. 1: Local density of states (LDOS) at the second vacuum 
layer above the 3d adatoms (2.86 A) on the Cu(lll) surface. 
The full lines refer to majority-spin states, the dashed lines 
to minority-spin ones. The two arrows show the protrusions 
discussed in the text: (a) corresponds to a split-off state, (b) 
is a d z 2 resonance. 



The calculated LDOS for 3d adatoms presented in 
Figffl exhibit several adatom-induced peaks. The LDOS 
refer to an unoccupied lattice position in the third vac- 
uum layer above the surface, i.e. 2.86 A above the single 
adatoms in the first vacuum layer. We focused on the 
region where the split-off state appears experimentally, 
i.e., around the bottom of the surface state, which in the 
actual calculations is located at —0.68 eV. We note that 
experimentally, the threshold energy is higher and is sit- 
uated at —0.45 eV—. This inconsistency is due to the 
LDA equilibrium lattice parameter we used since a test 
calculation with the experimental lattice parameter gives 
a value of —0.49 eV for the threshold energy. 

Let us start with a Cu adatom. Below E^, Fig^shows 
that two states appear in the LDOS. The first one (see 
arrow (a) in Fig^l is a split-off state situated at the bot- 
tom of the surface state (« —0.68 eV) as was already 
found by Limot et al^- and by Olsson et al£. To un- 
derstand the origin of the second protrusion (see arrow 
(b)) we plot in Fig|3the d-partial LDOS of the adatoms 
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which shows that it comes from a resonance of the d z 2 
state at 
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FIG. 2: d-contributions to the Local density of states (LDOS) 
of 3d adatoms on Cu(lll) surface. The full lines describe 
the majority-spin states and the dashed lines describe the 
minority-spin states. Since the Cu adatom is non-magnetic, 
the majority and minority virtual bond states coincide. 



For all the 3d-adatoms examined, we find for both spin 
directions a peak in the LDOS in the region below the 
threshold value of the surface state. However, we also 
find peaks at other energies, e.g. for Cu at —1.7 eV as 
mentioned earlier and for Co at —1.5 eV for spin- up and 
at around eV for spin-down states. In order to under- 
stand this behavior, we have in first place to understand 
the electronic and magnetic properties of the impurities 
and then which of the impurity states penetrate well into 
the vacuum, such that they show up in the LDOS of the 
second layer and can be detected by the STM. 

Let us first address the second question. In the vacuum 
region states with small in-plane fc-components fei | decay 
most slowly as a function of the perpendicular distance 
z from the surface. In fact, for a given energy E = —\ 2 
below the vacuum barrier, a wave function with in-plane 

component fen decays as e (x W' z . Therefore states 
with small fc||-values decay slowest; for higher fen part of 
the kinetic energy is "absorbed" by in-plane oscillations. 
Since the states with h i = show no in-plane oscillations 
and exhibit the full symmetry of the surface, we find that 
for the (111) surface states with s, p z and d z 2 show a 
slow decay in the vacuum region and can be well seen in 
STM, while other p- and d-states are strongly attenuated. 
An analogous argument holds for single adatoms on the 
(111) surface, since only these states exhibit the full point 
symmetry of the adatom-on-surface configuration 9 , and 
have thus no oscillations to absorb part of the kinetic 
energy. 

To understand the magnetic properties of the adatoms, 



we have plotted in Fig|2] the ^-contribution to the local 
density of states (LDOS) at the adatom site. For the Cu- 
adatom, we find a sharp d-peak at —1.7 eV, i.e. at the 
edge of the bulk d-band of Cu. This is a consequence of 
the repulsive potential, which the Cu adatom experiences 
in the first layer, shifting the d states to higher energies 
than in the bulk. All other impurities are magnetic and 
exhibit sizable moments, leading to a spin splitting of the 
so-called virtual bound states. For the different adatoms, 
the calculated local moments M s are: Cr (4.06 /is), Mn 
(4.28 fj, B ), Fe (3.21 (j, B ), Co (1.96 (i B ) and Ni (0.34 (j, B ). 
Note that the spin splitting is roughly given by / • M s , 
where / is the exchange integral of the order of 1 eV 11 . 
In the case of the Cr adatom, the minority peak is at 
higher energies and cannot be seen in Fig[2 while for Fe 
and Mn the majority peaks are at lower energies. 

Let us now come back to the interpretation of Fig^ 
showing the LDOS in the second vacuum layer, at the 
position above the adatom, for majority and minority 
electrons. Independently of the peak structure, we ob- 
serve a general increase of the DOS at higher energies, 
which arises from the increase of the spatial extent of the 
wave functions for larger energies. For Cu, the peak at 
— 1.7 eV coincides with the <i z 2-peak in the local DOS of 
the adatom, shown in Fig|21 For clarification we show in 
FigEl the local s-DOS and the d z 2-DOS of the adatoms 
(the latter reduced by a factor of 10). For Cu, as well 
as for the majority states of Co, we see a maximum and 
minimum in the s-LDOS at the <i z 2-peak position, aris- 
ing from the Fano-like resonant scattering of the s-states 
at the d z 2 -resonance. This effect cannot occur for a sin- 
gle adatom in free space or in jcllium, since the s- and 
d-orbitals are orthogonal. Therefore it is brought about 
by the reduced symmetry, i.e. the scattering at the sub- 
strate atoms. The LDOS-peak below —0.68 eV (arrow 
(a)) is the split-off state of the Cu adatom, induced by 
the attractive nature of the adatom potential in the first 
vacuum layer. The same states are also seen for the Co- 
adatom, more or less identical for both spin directions. In 
addition we see for the minority s-state of Co a Fano-like 
resonance behavior at the Fermi level, arising from the 
interaction with the minority d z 2 -virtual bound states. 

In the case of the Ni adatoms the virtual bound states 
for the two spin directions are only weakly split and more 
or less coincide with the energy level of the split-off sur- 
face state. Therefore in the local s-DOS both effects, the 
formation of the split-off state and the resonant scatter- 
ing at the d z 2 resonances, cannot be distinguished. How- 
ever, in the vacuum region (Fig^l the spin splitting of 
the majority and minority d z 2 -states can be clearly seen. 

For Fe-, Mn- and Cr- adatoms another effect can be 
seen in the vacuum LDOS (Fig^l and the local s-DOS 
in Fig|3J The intensity in the majority split-off surface 
state is considerably higher than for the minority state. 
This can have several reasons. For instance, due to the 
exchange splitting the majority potential is somewhat 
stronger than the minority one, leading to a smaller lat- 
eral extension of the split-off states and to a larger inten- 
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FIG. 3: (color online) Focus on the s-partial LDOS (black 
curves) of the impurity atoms and d z 2 -partial LDOS (red 
curves) reduced by a factor 10. Full lines represent major- 
ity spin, dashed lines minority spin. 




FIG. 4: (color online) Variation of the scattering length a 
versus the potential. Full lines correspond to a negative po- 
tential while blue one describe the case of positive potential 
(always repulsive). The figure is obtained with a small posi- 
tive energy value for an elastic s-scattering by a rectangular 
spherical potential depth. 



sity on the adatom site. In particular in the case of Cr, 
also the resonant interaction with the impurity d z 2 vir- 
tual bound state becomes important, strongly increasing 
the majority intensity on the impurity site as well as in 
vacuum. 

Moreover, we have noticed a small peak appearing in 
the Mn- and Fe-adatom majority d-LDOS (see arrow (c) 
in FigEJ) at the same position where the virtual bound 
states of Cu-adatom is situated (i.e. sa —1.7 eV). How- 
ever these protrusions have a different origin since there 
is no peak at —1.7 eV in the d z 2-LDOS for Mn- and 
Fe-adatom contrary to Cu- or Co-adatoms (see Fig|3J|. 
They appear at the remaining e?-partial LDOS (d xy , d yzi 
d xz and d x < i_ y i) which are strongly damped in the vac- 
uum. For symmetry reasons they do not hybridize with 
the s-LDOS explaining thus why we do not see a peak at 
the s-LDOS of Mn and Fe-adatoms contrary to Cu and 
Co-adatoms. We believe, however, these peaks can be 
interpreted as split-off states from a surface state at the 
M points, which shows a negative dispersion, such that 
a repulsive impurity potential leads for these d-states to 
a split-off state at higher energies, i.e. above the corre- 
sponding surface band. 

sp-impurities as adatoms: We consider now some 
sp-impurities as adatoms, for which the behavior is not 
complicated by the spin polarized e?-states. As first can- 
didate we consider Ca at the beginning of the 3d-series. 
The calculations give a well defined split-off state at the 
same position as for Cu and the 3d impurities, i.e. be- 
low the minimum of the surface band. In addition we 
perform calculations for Zn, Ga and Ge adatoms. The 
split-off state is still seen for Zn but not anymore for 
Ga and Ge atoms. The reason for this is that in the 



LDOS of Ga and Ge adatoms the s-states have moved 
to lower energies below the surface state minimum. In 
this case the s-scattering at the adatoms becomes effec- 
tively repulsive, so that no split-off state occurs. To 
explain this we note that the scattering behavior of a 
scattering center is directly related to the t-matrix, and 
only indirectly to the potential. Our results are in con- 
tradiction with the usual statement that any attractive 
potential leads to a split-off state of a two-dimensional 
surface state. This is not correct in our case, since the 
t-matrix of the adatom is basically a three-dimensional 
quantity. For s-scattering the t-matrix is for low energies 
E ~ given by the so-called scattering length a, being 
discussed in many books on quantum mechanics. The 
quantity a is the length where the extrapolation of the 
asymptotic form of the wave function for E — van- 
ishes. For the simple model of spherical potential well of 
depth V and radius ro the scattering length a is plotted 
in FigQ] For a repulsive potential the scattering length 
is positive and approaches the well radius ro for large V. 
For a weakly attractive potential a is negative. However 
when the potential V becomes stronger attractive, the 
scattering length assumes more stronger negative values, 
until it jumps. At a critical strength V = Vq from — oo 
to +oo, and is positive for further increased V- values. At 
the critical strength Vb a bound state appears at E = 0, 
moving to lower energies for further increased negative 
V values, and making the scattering length positive as 
for a repulsive potential. To compare with the real sit- 
uations of the adatoms, the potential of the transition 
metal atoms is sufficiently weak, that a split-off state ex- 
ist, since the atomic 4s-level is far above the Fermi level, 
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and the scattering length is negative. However for the 
Ge adatom the s-level has moved below the minimum of 
the surface state, so that the scattering length is posi- 
tive and no split-off state exists. When progressing in 
the atomic table a split-off state can only appear again, 
when in the next series of the elements, say from Rb to 
Ag, the 5s-level has moved down towards the Fermi level. 
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FIG. 5: (color online) LDOS in vacuum at the second layer 
above the impurities which are sitting in the surface layer. 
One can see the appearance of the split-off state above Ge 
(red line) impurities but not above Co (black line). Above the 
last one, a protrusion appears at w —1.3 eV on the majority 
spin channel which are due to the d state of the Co adatom. 



Impurities in the surface layer: The scattering of 
the surface states at impurities in the first surface layer 
is basically different from the scattering of s-electrons, 
since the effective potential for scattering is the difference 
between the potential of the impurity and the potential 
of the substituted Cu atom. Therefore all 3d-impuritics 
in the first layer do not show any split-off surface state, 
since the potential difference is very small (and moreover 
slightly repulsive). On the other hand for Ge impurities 
in the first layer, the calculations yield a split-off state, 
as is shown by the small peak in Fig(5J Apparently the 
difference in potential is sufficiently attractive, such that 
a weakly localized state is formed. Therefore we obtain 
the opposite trend as for the adatoms. Transition metal 
impurities exhibit a split-off surface state as adatoms, but 
not as substitutional impurities in the first layer, whereas 
for Ga and Ge just the opposite is true. 



many cases below the Kondo temperature; e.g., a char- 
acteristic Kondo feature in the spectra was observed for 
Co on Cu in Re& 

The Kondo effect is characterized by a narrow 
Abrikosov-Suhl resonance of the DOS at Ep which is 
absent in our calculations. However, the split-off states 
are well below Ep. Furthermore, the Abrikosov-Suhl ap- 
pears by the interaction of impurity s-states with the 
surface band. Therefore our results on the split-off state 
are physically relevant. 

On the other hand, the spin-dependent spectra of the 
magnetic impurities should be corrected towards an av- 
eraging of the two spin channels, if the temperature is 
below the Kondo temperature. Although the Kondo fluc- 
tuations kill the magnetic moment, the splitting of the 
d virtual bound states remains, corresponding to singlc- 
and double- occupancy of the impurity local state. 



V. SUMMARY 

We have performed self-consistent calculations on sin- 
gle impurities deposited on Cu(lll) surface in order to in- 
vestigate the split-off state recently seen in experimenljiii. 
We show the existence of two kinds of state localizations: 
One is due to an attractive potential below the V surface 
state and the other one is due to a repulsive potential 
above the M surface state. We found that Ca, all 3c? and 
Zn adatoms produce the first resonance which is not the 
case for the sp adatoms Ga and Ge even if fundamentally 
it is known that any attractive potential should lead to 
a split-off state at the bottom of a two-dimensional elec- 
tron state. Its presence for Ca means that s states have 
a stronger contribution to its realization then d states. 
The behavior is totally different if the impurities are em- 
bedded in the surface layer. In particular, Ga and Ge 
lead to a protrusion at the bottom of the surface state. 
One question addressed in this work is how to explain 
the non-existence of state localizations for sp-adatoms 
and 3<i-impurities embedded in the first surface layer. In 
the latter case, the reference potential is then the Cu 
one which is not too different from the potentials of all 
3c? impurities. Therefore we do not see any localization. 
The case of sp adatoms is explained by the sign change 
of the scattering length for attractive potentials, which 
then may act effectively as repulsive potentials. 



IV. LIMITATIONS OF THE LSDA 

We no discuss the limitations of our calculations, in 
particular concerning the Kondo effect which cannot be 
captured by the LSDA. It is well-known that, at low tem- 
peratures, the spin moment of the magnetic impurities 
fluctuates, so that these appear non-magnetic. The tem- 
peratures at which the experiments are conducted are in 
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